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Abstract 



We demonstrate a practical and efficient metliod for generating toric Calabi-Yau quiver the- 
ories, applicable to both D3 and M2 brane world-volume physics. A new analytic method is 
presented at low order parametres and an algorithm for the general case is developed which has 
polynomial complexity in the number of edges in the quiver. Using this algorithm, carefully im- 
plemented, we classify the quiver diagram and assign possible superpotentials for various small 
values of the number of edges and nodes. We examine some preliminary statistics on this space 
of toric quiver theories. 
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1 Introduction 

The works of [1,2] have brought about tremendous interest in the M2-brane, and following [3], where 
the world- volume theory was realised as a Chern-Simons quiver theory in (2 + l)-dimensions, much 
activity ensued (q.v. e.g. [4-23]). Even though nice parallels are drawn between the familiar case 
of D3-branes probing Calabi-Yau threefold singularities and the present circumstance of M2-branes 
probing Calabi-Yau fourfold singularities, the latter situation is far less understood. Complications 
arise in the various analogues of the wealth of properties enjoyed by the D3-brane, such as singularity- 
resolution in relation to (un)Higgsing, the "inverse algorithm" for systematically constructing the 
world-volume gauge theory, Seiberg duality, etc. Part of the issue arises from the new possibility of 
turning on G-fluxes on torsion-cycles in the dual AdS geometry [20,21,23]. 

Confronted with these seemingly untamable complexities, an optimistic path, guided by brane 
configurations, had been trodden. It is by now well-established that the (3 -I- l)-dimensional su- 
persymmetric quiver gauge theory of D3-brane probing a toric Calabi-Yau singularity is completely 
described by a planar, periodic tiling of NS5-branes and D5-branes, or, equivalently, by an auxiliary 
bi-partite graph on a torus known as a dimer model. It is suspected that this persists to the case of 
our (2 + l)-dimensional quiver Chern-Simons theories, a conjecture fortified by the observation that 
almost all such theories discovered to date, being of the M2-brane world- volume, descend from parent 
D3-brane theories. More precisely, take the quiver and superpotential of a world-volume theory aris- 
ing from a D3-brane probing a Calabi-Yau threefold, add Chern-Simons levels and terms respectively 
to the graph labeling and the interaction, the new moduli space is computable to be a Calabi-Yau 
fourfold, apparently rather different, but is expected to be the space which a M2-brane probes. This 
"forward algorithm" is succinctly presented in [17]. 

Led by this beacon, a systematic, taxonomic study of toric quivers was undertaken in [17] and 
it is in this spirit that we wish to proceed with our current investigations. Our purpose is twofold. 
First, we wish to continue with [17] and exhaustively and progressively classify toric quiver diagrams 
with superpotential, whereby establishing, for each given number of nodes and edges, all possible 
parent (3 -I- l)-dimensional theories which could descend to (2 -I- l)-dimensional Chern-Simon quiver 
theories. We will see this is facilitated by a tremendous improvement of generating quivers, using 
an efficient notation and algorithm. Second, we would like to take advantange of this new tool and 
attempt an exploration of the space of quiver theories along the vein of [34,35], and gather data for 
our experimental probe of the string geometric landscape. 

The paper is organised as follows. In the next section we briefly review some topics in graph 
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theory, complexity and toric Calabi-Yau quiver gauge theories. This is followed by a restatement of 
the problem, rephraseded in this language. Section 3 offers a new method by which toric Calabi-Yau 
quiver theories can be generated, a new notation for specifying quiver diagrams and an example of 
both of these in the form of the complete generation of a low order set of quiver diagrams. Section 
4 presents an algorithm used for generating complete sets of quiver theories. In section 5 we show 
our results, presenting some complete sets of quiver diagrams and associated superpotentials; as well 
as some statistics of the total numbers inequivalent quiver theories for some low order parametres. 
Finally, in section 6 we conclude with a discussion of the work, and its futher prospects. 

We would like to draw the reader's attention to a beautiful parallel work by [36] , which also classifies 
toric quiver theories, but from the point of view of brane tilings. The two works thus complement 
each other harmoniously. 

2 The Taxonomical Problem 

We start with introducing the concept of a toric Calabi-Yau quiver theory (TCYQT) in an abstract 
sense. We refer the reader to [29], which first brought the study of quiver theories to gauge theories, 
for details, as well as to [31] for a rudimentary review. 

2.1 Some Rudiments on Graphs and Complexity 

First, let us recall some elementary facts on graph theory and algorithmic complexity, which will be 
crucial to our constructions later. A graph is an ordered pair {G,E), where G is a set of nodes and 
E is a set of pairs of nodes, which represent the two ends of an edge. In a directed graph, E contains 
ordered pairs, which tell us the edge direction [24]. Node i is adjacent to node j if the directed edge 
is in the graph. 

A path on a directed graph with nodes G = {ni,i = {l..p}} is an ordered set of edges {E = 
{ej,j = {l..q}} where ej — [nj^^rij) and Cj^i = (rij.ni). That is, we follow the arrows on the graph 
between two given nodes k ^ I. In a directed graph, each node has an 'indegree^ and 'outdegree\ 
which are the number of edges entering and leaving the node respectively. A connected graph has a 
path between any two nodes and a cycle is a path which ends at its starting node. A graph is Eulerian 
if there exists at least one cycle which visits every edge once and only once [24] . A quiver graph can 
be viewed as an extension of a directed graph and is the quadruple {G, E, d, R) where d is a set of 
node labels which tell us the gauge groups of the theory and R is a. set of algebraic labels for the 
edges. The choice of d and R is known as the quiver representation [25]. 

A major motivation of this investigation is to establish an efficient method of explicitly constructing 
all toric quivers. The best way to analyse the quality of an algorithmic method is to examine its 
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complexity, i.e., how the time the algorithm takes to run increases with the size of its arguments. An 

algorithm has polynomial complexity [27] if its running time increases as some power of the size of the 
argument. For example, looping through a list has complexity 0{N), where N is the muiibcr of terms 
in the list. As a general rule, algorithms with polynomial complexity are known as "fast" algorithms 
and ones with exponential complexity arc "slow" [27]. 

There is an important area whore graph theory and computer science meet and that is the question 
of graph isomorphism. Two graphs are said to be isomorphic, and so belong to the same graph 
isomorphism class if there exists a relabeling of the nodes which makes the graphs identical. Currently 
there is no known simple property of a graph which can be used to check whether it is isomorphic 
to any other. It is an open question as to how efficient an algorithmic solution for checking graph 
isomorphism can be. The current quickest solution to the isomorphism problem is found in [28]. For 
a detailed discussion of the problem see [26]. 

2.2 Toric Calabi-Yau Quiver Gauge Theories 

Prepared with the above elements, we now introduce quiver graphs in the context of brane world- 
volume gauge theories. A quiver theory is a pair (Q, W) where Q is a finite quiver diagram (directed, 
labeled and represented in the aforementioned sense) and is a formal polynomial in the arrows of 
Q, called the superpotential. The finiteness signifies that there is a finite number of nodes, which 
are to correspond to gauge factors in a gauge theory, and a finite number of directed arrows, which 
are to correspond to bi-fundamental fields. The labelings are integers to be assigned to the nodes, 
and are to denote the rank of the gauge factors as well as possible Chern-Simons levels. Let Q have 
G nodes, E arrows and let W have Nt terms. The arrows will be denoted as X^/^, signifying fields 
charged under nodes a and b, with possible multiple arrows indexed by i. The superpotential is then a 
polynomial in these fields, each monomial term is a gauge invariant and is thus a product over closed 
loops in the quiver; furthermore, we require the monomials to be at least cubic since quadratic terms 
are mass terms which could be integrated out. For now, we shall take all fields to be simply 
complex numbers, this will be the first of so-called toric conditions, and in particular corresponds to 
the case of a single brane-probe. In general, we should promote these to matrices, in which case there 
are multiple parallel coincident branes. 

We shall represent Q by the incidence matrix d, which is a G x matrix of ±1 and entries: 

each row enumerates a node and each column, an arrow, such that an arrow from node a to b will 
have, in the column representing this arrow, the a-th row being —1 and the 6-th row, all other 
entries are 0. Clearly, the sum over all the rows gives a row of zeros by definition. Note that this 
matrix does not capture so-called adjoint fields which go from a node to itself. 

To each such a theory {Q,W) one can construct a so-called vacuum moduli space (VMS). This, 



4 



by construction, is the geometry which a branc should probe. The procedure by which one finds the 
VMS from the pair (Q, W) has been called the "Forward Algorithm" [32]. In the particular case where 
the VMS is a toric Calabi-Yau space, this algorithm can be succinctly summarised as manipulation 
of integer matrices and combinatorics (sec the flow-charts on page 7 of [16]). The reverse problem, 
of geometrically engineering the brane world-volume theory given the vacuum geometry is likewise 
called the "Inverse Algorithm" . 

Having introduced the above fundamentals of the quiver theories which we will systematically 
study in this paper, we now move on to placing some mild constraints so as to restrict ourselves to a 
sub-class which had been of primary interest over the last decade, viz., toric Calabi-Yau conditions. 
As we mentioned above, we shall first take all gauge groups to be U{1), whence all fields are complex 
valued. Next, we shall require that the VMS, computed from the Forward Algorithm, be Calabi-Yau, 
this requires that the sum over the columns also gives a column of zeros. Now, in (3 + l)-dimensions, 
this is a consequence of anomaly-cancelation: that for each node, the indegree and outdegree should 
equal. Note that this is a sufficient but not necessary condition for the VMS to be Calabi-Yau, 
however, for simplicity we will henceforth impose it on all our quivers. 

In order that the VMS be toric Calabi-Yau, further constraints are to be met by the superpotential. 
First, as pointed out in [32] and dubbed the "toric condition" , each field must appear exactly twice 
and with opposite sign. Second, it is now well-established that toric Calabi-Yau brane theories admit 
a tiling or dimer description (see [33] for excellent reviews); in order to admit a tiling, we must also 
adhere to a topological condition 

G-E + Nt = 0. (2.1) 

2.3 Summary of Constraints 

We are thus confronted by a classification problem. Our strategy will be to generate incidence matrices 
of size G X E which represent the quiver diagrams, such that the points 1 and 2 of the following set of 
conditions are met. To each such a quiver Q we must generate all possible superpotentials W which 
satisfy conditions 3 to 6: 

1. Every column in the incidence matrix must sum to zero (incidence condition). 

2. Every row in the incidence matrix must sum to zero (Calabi-Yau condition). 

3. Every superpotential must satisfy the toric condition. 

4. Every superpotential must contain x = E — G terms (tiling condition). 

5. Every term in the superpotential must correspond to a loop in the quiver diagram (gauge 
invariance). 
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6. Every term in the superpotential must contain at least three fields. 

Computationally, we may proceed as follows. The initial space of possible matrices are those 
which contain entries that are +1, —1 and 0. Conditions 1 and 2 should be imposed on them, 
producing a quiver diagram incidence matrix and reducing the space of possible incidence matrices. 
Then a representative of each graph isomorphism class is required from this reduced set. The brute 
force approach to generating possible incidence matrices has at worst complexity 0(3^'^^^^-'), as each 
element in the matrix can be one of three entries. Next, we also wish to find a constructive way 
of forming supcrpotcntials. The strategy will be to find a way of producing a superpotential using 
the given quiver as a foundation. Here, the brute force algorithm for superpotential generation has 
complexity 0(2^ ), as the set of all sets of monomial terms is required. Clearly, the brute-force method 
is not feasible and from a computational perspective a fast algorithm is desireable, though due to the 
difficulty of the problem the emphasis is on finding a simple workable solution, easily testable to 
eliminate errors. We will now devise such an algorithm. 

3 New Notations and Generations for Toric Quivers 

As we saw above, we need an efficient, or at the least, a practical method of systematically generating 
toric quiver theories. Preliminary attempts were made, using standard partition routines, especially 
implemented on Mathematica® , in [16]. However, there, the taxonomical study was already hindered 
by the growth rate of the complexity. In this section, we will start afresh, and introduce a new notation 
for generating quivers, engendering all from a polygonal shape, and then rather effectively enumerate 
our desired theories. 

3.1 Collapsing Polygons and Constructing Super potentials 

We begin with the simple observation that every toric quiver diagram admits at least one Eulerian 
cycle since every node has equal indegree and outdegree due to the Calabi-Yau condition. The original 
proof of this is due to Fleury [24] and a proof is included in Appendix 1X1 This realisation that a quiver 
diagram admits an Eulerian cycle gives a route to a smaller space of possible toric incidence matrices. 
Consider a polygon, with nodes at the corners and clockwise directed edges. By pinching together 
two nodes in the polygon, it is collapsed to a graph with fewer nodes. This graph clearly obeys all 
conditions for a toric Calabi-Yau quiver diagram. 

This above intuition can be formalised. Suppose we wish to study a toric quiver with G nodes 
and E arrows, we define a; = i? — G as an order parametre as was done in [16]; this of course, will 
turn out to be the number of terms Nt in the superpotential by condition (|2.ip . Now, we can define a 
graph operation C{i,j) on the polygonal quiver as one which pinches the graph between nodes i and 
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j. After this is applied, the resulting node is labeled simultaneously with both i and j. The polygonal 
quiver together with this pinching operation is shown diagrammatically in Figure [T] 



Q o 




Figure 1: Pinching a cyclically directed polygonal graph by the operation C{i,j). 

It turns out that the above simple concepts already suffice to give us a powerful result. 

PROPOSITION 1. By drawing an E-sided directed polygon with G numbered vertices and operating 
with C{i,j) X = E — G times, for different pairs {i,j), every unique directed toric quiver diagram is 
generated. 

Proof. : Let the set of toric quiver diagrams with E — G ~ x he Qu^ , the set of directed polygon 
graphs be P and the set of graphs formed by acting on P with C{i,j) x times be denoted CP^ . We 
now prove that Qu^ equals CP^. In one direction, consider an i?-sided directed polygon. If C{i,j) is 
applied to nodes i and j the resultant node has two edges entering and two edges leaving the combined 
node which we denote as («, j). Repeated application of C((«, j), /c) on the resultant node adds another 
edge entering and leaving each time. So the graph still satisfies the Calabi-Yau condition. Applying 
G{i, j) X times with different (i, j) constructs a directed graph which satisfies the Calabi-Yau condition 
and still prodces toric quiver diagrams, with E edges and E — x = G nodes. Therefore CP^ C Qu^. 

Next, consider a quiver, and walk an Eulerian cycle through it. Every time we reach a node with 
multiple edges entering and leaving we split the node. The edges we travel along are now the only 
edges attached to the new node. This is the inverse operation to C(i, j). This process forms a directed 
graph where E = G, and constitutes a member of P. Thus, Qu^ C GP^ . Therefore, 

CP^ C gu^ & Qu^ C CP"^ ^ Qu^" = CP'' 

□ 

A map TT : P i~> Qu^ can now be defined as a way of producing toric quivers. This is a way of 
guaranteeing a very large reduction in the space of possible incidence matrices. 

Having constructed the toric quiver Q, we now assign the superpotential W. 
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PROPOSITION 2. Toric superpotentials can be constructed by decomposing the quiver diagram 
into disjoint loops to generate terms and then matching up the collections of terms so that the total 
number of terms is x. 

Proof. Drawing all loops which correspond to terms of one sign on top of the quiver diagram creates 
a collection of loops which contains every edge once. Thus every set of disjoint loops which covers the 
quiver diagram represents a collection of terms. We can then take two collections, one with n terms 
and one with x ^ n and subtract one from the other. On the condition that the two have no common 
terms we have made a superpotential with x terms. Doing this for all possible combinations of term 
collections will generate all toric superpotentials. □ 

The term collections are constructed by finding all ways of splitting the quiver diagram up into 
disjoint cycles. Consider walking through the quiver, there is a choice of where to go only at the 
contracted nodes. We define the concept of edge assignments at each node by a set of rules EAi = 
{{Ea,Eb) : Ea enters i and Et leaves i}. By constructing all EAi V i: indegree(i) > 1 we produce all 
sets of disjoint loops of the quiver diagram. 

To illustrate, let us use a concrete example. Let us study the (G, E) — (3, 6) quiver diagram. Note 
that we will henceforth use this pair notation to denote quivers. We have drawn the (3, 6) quiver in 
Figure [2] and now will produce all toric superpotentials for it. Of course, this is only a demonstrative 
example, since there will be2: = 6 — 3 = 3 terms in the superpotential and we will see there are 
non-minimal cycles. In the actual physics, we may need to promote the node ranks to non-Abelian 
and shuffle traces amongst the operators. All sets of two loops that can be used to form the gauge 



Figure 2: For the {G,E) = (3,6) quiver, we construct all two term collections. One term is the product 
of the fields in the dashed loop, the other is the product of fields in the solid loop. These will then appear 
in the superpotential, with opposite signs. 
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invariant terms are shown in Figure [51 As well as these loops, we also have the term made from an 
Eulerian cycle of the quiver, namely Taii = X2'^ X^]^ X^^\ looped over all the arrows. 

The term collections are then matched up to give x = 3 terms and hence construct all superpotentials: 

Wi — A;^2 ^23 ^31 +^12-^23-^31 - -l^ all 

M/ - y(2)T^(l)T^(2) rp 

W2 — A;^2 -^23 ^31 +^12-^23-^31 ~ all 

W - 1^(2)^^(2) „ 

— -^12-^23-^31 +-^12-^23-^31 - ^ all 

W - y(2)-^(2)„(l) „(l)y(l)y(2) „ 

W4 — A;^2 -^23 -^31 +-^12-^23-^31 —-''all- [o-^) 

Indeed, the above combinations are exhaustive and we will perform similar such partitions for 
all the toric quivers which will be constructed later. Now, a limit on the quivers which admit toric 
superpotentials can also be found. The quiver must be able to be decomposed into at least ceil(|) 
loops (we use the standard notation that ceil(x) rounds x up to the nearest integer and floor(x) rounds 
X down), as this is the minimum number of terms required from one term collection to give x terms. 
Given that there must be at least three fields in each term, the maximum number of terms is E/3. 
Thus, 

E > 3ceil(|) . (3.2) 

3.2 AB notation 

Now, we introduced the map tt : P ^ Qu^ above which explicitly constructs toric quiver diagrams 
from polygonal contractions. To specify the maps we propose a notation which depends as little as 
possible on labeling the nodes and which reflects the symmetry properties of the quiver. This will 
turn out to be important for redundancy elimination. 

We take a walk through a member of P to build the quiver and make the following definitions: 

• Let C{a) :— C{i,i + a) where i is the node at which we are currently stopping and a is the 
number of steps ahead. 

• Let W{b) be the operation of moving along b edges through the nodes in their labeled order 
from the current position. 

• We call the process C{a) which contracts two nodes which have already been contracted as a 
redundant contraction. 

• We call a node which has been contracted with another as a c-node. 

We can now define our walk as follows. Due to the nodes being all labeled by 1 in the toric Abelian 
case, the walk can begin at any node that will be contracted. First, apply C{ai) and then W{bi), 
where ai is the number of edges in the first contraction and bi is the number of edges to the next 
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node we wish to contract. Then, apply C(a2), followed by 1^(62), and so on. Stop once x contractions 
have been performed. Finally apply W{hx) to return to the starting node. Thus any quiver diagram 
can be represented by two ordered sets of numbers: 

A = {fli : i = 1, ..,x] 

B = {b,:i = l.-,x} , (3.3) 

where the a,; in the A define the contractions in the order in which they are applied, and the bi in the 
B define the walks between contractions. Therefore, a quiver can be described by applying elements 
from the A and B in turn. That is, using ai followed by hi then 02 and so on. Thus a general AB is 
a constructive map: 

(A, B):P^Qu'. 




Figure 3: The action of the constructive AB maps on the polygonal directed quiver to a toric quiver. 
The red node is the current location of the walker, the left hand column represents the application of the 
contractions {ai), from the red node to the green node. The right hand column shows the state after the 
application of the walks {hi). The diagram shows the application of A = {2, 2}, B — {1, 4}. 

This AB notation has a number of useful symmetry properties which are due to A and B containing 
relative node positions on the quiver. Rotating the quiver amounts to a cyclic permutation of the 
elements of both A and B and reflecting the quiver diagram is also manifest, but in a more complicated 
way. We refer the reader to Appendix |X] for more details on these symmetries. It is important to 
note that the notation is much more efScient than an incidence matrix, requiring only 2x numbers to 
label a quiver. Using these concepts and this notation we can restrict the space of possible incidence 
matrices and produce an algorithm much more efficient than simple brute force. 
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3.3 Warm-up: The x = 2 Problem 

As an example for applying these above methods, let us classify, as a first example, the x — 2 quiver 
diagrams, that is, all quiver diagrams formed by two contractions from some polygon. These are the 
two-term quiver theories studied in [16] and are quite simple quivers with A — {ai, 02}, B — {6, E~b}. 
With the above ideas, it should be possible to analyse this case easily. For simplicity only quivers with 
no self-adjoint loops are considered, as x — 2 quivers with self-adjoint loops are trivial extensions of 
X — and x ~ 1 quivers, i.e., we just add a loop to a single node for the a; = 1 quivers, or two loops 
for the a; = 0. 

Each X = 2 quiver is made by applying some contraction operation to an a; = 1 quiver, so we must 
consider the possible x — 1 quiver diagrams first. We need to work out which quivers are unique. One 
immediately sees that applying C(a), where a is the node spacing, will give us an a; = 1 quiver, however 
if a > fioor(i?/2) we can flip this diagram to become a quiver with a' = a — fioor(i?/2). Therefore 
all the X — 1 quivers are formed by applying the operation C{ai) for ai — 2,3,..., floor(£'/2). To 
produce the x — 2 quivers we first apply the walk operation W{bi), to move to another node in the 
quiver. Then another contraction C{a2) is applied . The values of bi and 02 which can be applied to 
produce non-isomorphic quivers split into a number of regimes, based upon their values. 



Figure 4: The solution to the warm-up problem of x = 2, by splitting up into regions 1 to 4. 

We present the procedure graphically in Figure 21 where we have labelled 4 regions. The idea for 
solving is that we will systematically contract all nodes with each other, but single out those which 
will produce isomorphic quiver diagrams. We first look within the first quarter of the general a; = 1 
diagram above, confining contractions to this quarter. We then match up all the contractions across 
quarters. As the quarters cut the diagram along potential symmetry lines, we can use this view of 
the diagram to easily work out which contractions will be isomorphic to one another. Due to the 
symmetry conditions, if we perform C(2, 6) on the above diagram, we have also performed C(3, 7). 
So we must factor this condition into our systematic categorisation of the possibilities. 

The matchings we choose to do go as follows: In region 1 we match up with regions 1 and 2. In 
region 2 we match with region 3 (this is isomorphic to matching region 1 with 3 or 4, and region 2 
with 4). In region 3 we match with regions 3 and 4. 

By working through these cases we will discover all the possible values of AB which let us create 
quiver diagrams not isomorphic to one another. We also note the no self-adjoint loop condition, this 
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has a number of effects. First, we restrict the members of A to bo strictly greater than one. Next, 
we remember in computing the conditions below that we do not want to create a self-adjoint loop 
by pulling any node which is one edge away from another together. Specifically, this can happen for 
a; > 1 quivers as applying contractions when there is a loop already in the diagram can lead to a 
situation where the a of the contraction is greater than one, but the distance of the node from the 
multiply connected node is only one. For the case below this is easy to avoid. However, more complex 
cases would prove a very real problem. 

First we have from the a; = 1 result that 2 < ai < floor{E/2), and let bm = ci + floor{{E — a)/2) 
to simplify notation, then let < 6 < 6m, and so: 

if ai ^ E/2 then we find following conditions on a2 for a given b and ai: 

if 6 = then 2 < a2 < floor{ai/2) 
if < 6 < floor{ai/2) then 2 < 02 < ai - 26 
if floor{ai/2) <b<ai then ai — b + 1 < a2 < bm — b 
if 6 = ai then 2 < a2 < 6m — 6 
if 6 > ai then 6m - 6 < a2 < - 26 + a ; (3.4) 

if ai = E/2 then we gain another line of symmetry and the only matchings required are: region 1 
with regions 1 and 2; region 2 with region 3. 

if 6 = : 2 < a2 < floor{ai/2) 
if < 6 < floor{ai/2) : 2 < 02 < ai - 26 
if floor{ai/2) <b<ai :a — 6+l<02<6m+a — 6. (3.5) 

Finally, we consider that we could have constructed isomorphic quivers which have the same loops 
but in the opposite order, i.e., with ai and 02 swapped. In order to make sure this does not occur, 
we restrict the value of 02, so that the loops always have a definite order depending on their size. 

a2 < ai 

a2 > 1 . (3.6) 

This exhaustively produces all possible quiver diagrams for .t = 2. 

It would be incredibly complex to conduct the same kind of analysis for larger x. Due to this 
difficulty in implementing this kind of constructive solution for arbitrary x, the next stop is to try an 
algorithmic solution, and use a graph isomorphism test to find quiver diagrams representative of their 
isomorphism class. To this we now turn. 
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4 The Algorithm for Generating Toric Quiver Theories 



Our warm-up example clearly lends itself to an algorithmic approach. In this section, we briefly 
outline this algorithmic procedure which we device to systematically generate the toric quivers and all 
possible associated superpotentials. The method has been implemented in C and is presently being 
organised into a convenient software. This section, together with the companion appenix, are rather 
technical and the readers interested only in the results can safely skip the ensuing text and move onto 
the next section where the results of our classification will be presented. 




Figure 5: Flowchart of the basic structure, showing the main processing sections and decisions of our algorithm 
for generating toric quiver theories. 



Our methodology is as follows. We will generate all the toric quiver diagrams by finding all possible 
AB, and test them to pick a representative of each graph isomorphism class. As AB is represented by 
a list of numbers we treat the problem as though the AB are numbers in base E that can be counted 
through; digits greater than E need not be considered as we would just loop around the whole diagram. 
Next, given a quiver diagram we wish to check if it is isomorphic to any quiver diagram which can 
be written as a "smaller" AB, that is, a smaller base E number. The 'smallest' AB is then kept as 
the representative of the class. With this quiver diagram we find the possible superpotentials using 
the method detailed in Section [3l The algorithm terminates when the AB, (interpreted as a base E 
number), overflows. The flowchart Figure [5] sets this out in detail. 
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4.1 Counting 



When implementing, the isomorphism test should check as few quivers as possible to maximise ef- 
ficency. The AB can be constrained to reduce the number that must be checked. To do this, the 
so-called ^moduli rules^ are used. These rules constrain the AB so the operation it represents never 
crosses the starting node and are given by imposing 



Moduli rules are used instead of restricting using symmetry since using symmetry leads to problems 
with redundancy and is more computationally complex. We can also place limits on the members of 
A and B. We require all a.; > 1 since C(0) is an identity and that 6^ > 1 to eliminate redundancy; 
again see Appendix |A] for proofs of the above statements. With these rules it is a simple matter to 
write the code which increments the AB, treating it as a base E number, and obeys these conditions. 
To implement the counting part of the algorithm a number of routines are used; these are detailed in 
the Counting section of Appendix IB. 21 

4.2 Isomorphism 

A simple graph isomorphism test would require the current quiver diagram to be checked against all 
quiver diagrams found so far, so a method that can test a single possible quiver at a time is preferable. 
A solution is to construct all possible ways of writing the quiver diagram in AB notation. This is 
done by walking all Eulerian cycles and counting the edges between c-nodes in these cycles to find all 
the AB. It is basically the reverse process of collapsing the polygon with the AB; we refer the reader 
to appendix [b1 Walks need only start from c-nodes, since C{a) is always applied first in an AB, 
increasing the efficiency. In fact, in Appendix [Bl we show that the routine has complexity 0{E'^^). 
The "smallest" way of writing the quiver diagram is then taken as the representative. We detail the 
routines in the Isomorphism check section of Appendix IB. 21 

4.3 Superpotential Generation 

In order to produce the superpotentials we first identify all c-nodes in the quiver diagram and then 
generate all edge assigments. For each edge assignment, we trace through the loops which it forms 
in the quiver, and then generate the gauge invariant terms. This is precisely the method discussed in 
section|3l As the order of fields within a term is irrelevant they are ordered in a consistent way to check 
if term collections share terms. Then the term collections are matched up to form superpotentials 
with the correct number of terms. The sub-routines are explained in the Superpotential Generation 
Section of Appendix IB. 21 
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5 Classification Results 

In this section, we put the above theoretical algorithms to practice and present the complete quivers 
for some low values of the order parametres. Many quivers have already eluded current methods of 
attack due to their complexity and our algorithm can generate them efficiently. Moreover, we will 
initiate an exploration of this space, consisting of many quivers new to the literature, and see what 
structures we can unravel. 

5.1 Quiver Diagrams for Low Order Parameters 

We now present some of the low order quiver diagrams for x = 0, 1, 2, 3. The a; = and x — 1 quiver 
diagrams are trivially an E sided polygon and single contractions A — a, B ~ E, so we need not draw 
them explicitly. Moreover, we will for now focus on what we call base quiver diagrams; that is, 
diagrams without self-adjoint loops. As we mentioned above, these adjoints are not captured by the 
incidence matrix because they would be ambiguously represented by a column of zero. One can be 
put these in by analysing candidates for the superpotentials as was done in [16]. 

The (G, E) equaling to (3, 5), (4, 6) and (5, 7) quiver diagrams are the lowest order x = 2 diagrams 
which are not trivially found. We see from inequality (j3.2p that if a: = 2, then E > 3. Requiring for 
now that there are diagrams with no self-adjoint loops takes this limit up to -E = 5. The (3, 5) and 
(4,6) diagrams match those analysed by [16]. 



Figure 6: {G,E) ~ (3,5) Quiver Diagrams, self adjoint loops are implied in the graphs with edges fewer 
than E = 5. 

Diagrams with self adjoint loops are not shown below (except for the (3, 5) case) as they are the 
trivial x = 1 and x = diagrams with one (or more) self adjoint loop(s) attached to some node(s). 
There are three (3, 5) base quiver diagrams shown in Figure [6] a single quiver diagram with no self- 
adjoint loops, one with a single self-adjoint loop and one with two self-adjoint loops. There are five 
(4,6) base quiver diagrams shown in Figure [7] (allowing for adjoints there would be 7). In Figure [8] 
the seven (5, 7) base quiver diagrams are presented (allowing adjoints there are ten (5, 7) quivers). 

Next, let us move on to a; = 3, and present some low order complete sets of quiver diagrams. This 
has not appeared in the literature before. We find that E > 6 from inequality (|3.2p and so present the 
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(3) -@ 



Figure 7: {G,E) ~ (4,6) base Quiver Diagrams, with no self adjoint loops. The diagram with red nodes 
does not admit a toric superpotential even if fields are not constrained to only take complex values. 




Figure 8: All {G,E) = (5,7) quiver diagrams with no self-adjoint loops. The diagrams with red nodes 
have vanishing toric superpotentials even if the fields are not constrained to only take complex values. 

three (3, 6) base quiver diagrams with no self adjoint loops in Figure [9l The quivers with sclf-adjoint 
loops are those in Figure [SI but with loops attached, giving six quiver diagrams in total. In FigurdTUl 
the five (4, 7) quiver diagrams with no self-adjoint loops are presented. The single self adjoint loop 
base quiver diagrams have already been shown in Figure [T] 

To demonstrate that the algorithm is still efficient for higher values, let us move on to an example 
set for x = 4, with (G, E) — (4, 8), which again is new to the literature. We present all base quivers in 
Figure [TT] and there is a total of 16. From this list we recognise, among others, the 8-arrowed quiver 

16 




Figure 9: {G,E) = (3,6) Quiver Diagrams, with no self adjoint loops. The diagram with green loops has 
a single toric superpotential and the diagram with red nodes admits no superpotentials. 



for the cone over the zeroth Hirzebruch surface. 
5.2 Assigning the Superpotential 

The above are simply quiver diagrams, providing any of them with a quiver representation gives a 
new physical theory with E matter superfields and E — x gauge groups. These classifications and the 
many others made by the algorithm for higher E and x therefore produce many new physical theories, 
the majority of which have not yet been studied and were previously unknown. It would certainly 
be interesting to study their VMS to identify the geometry. Some, of course, have appeared in the 
literature in interesting and significant applications. The single (2, 2) quiver diagram produced by the 
program was first used as for the supersymmetric gauge theory associated to the conifold [30] . One of 
the (3, 6) quiver diagrams was one of the very early quiver theories, used as a model for a conformal 
supersymmetric gauge theory [29]. These diagrams, among others, confirm the results found using 
the program and show that so far we have barely scratched the surface of the space of toric quiver 
theories. 

To produce full TCYQTs we require toric superpotentials. The algorithm has found superpoten- 






Figure 10: {G,E) — (4,7) Quiver Diagrams, with no self adjoint loops. 
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Figure 11: {G,E) — (4,8) Quiver Diagrams, with no self adjoint loops. 

tials for a large number of quiver diagrams. We again systematically work through the low order 
diagrams. The x — 2 quivers require two terms in the superpotential. In the TCYQT case both the 
positive and negative terms contain all fields and so will cancel giving no toric superpotentials in the 
Abelian case of a single brane for any x — 2 quiver diagram. The superpotentials for the (3, 6) quivers 
are found to be identical to the example results in[31 as they should be. It is also found that the (3, 6) 
quiver diagram with green nodes in Figure [9] admits a single superpotential: 

W = Xj^2-''^23^31 + -''^1 2-^^23 ^31 ^ -'^12 -'^23 ^3 1^1 2-^^23 ^31 ■ (^-l) 

The (3, 6) diagram with red nodes does not admit any superpotentials which do not vanish for the 
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Abelian case. Including the quivers with self-adjoint loops there are twelve (3, 6) TCYQTs in total. 




Figure 12: An example of {G,E) = (4,7) quiver diagram, for which we present all superpotentials in 
Equation ((O)) . 

As a more detailed example we now stdy a single (4, 7) quiver diagram, presented in Figure 1121 
which has no self-adjoint loops and admits multiple toric superpotentials: 

W - ^(1)^(2)^^(1) -^(2)^(1) ^^(1)^^(1) „ 
14^1 — A;^2 -^23 ^31 +-^12-^23-^34-^41 ^ -l^ all 

W - ^(1)^(1)^(1) -L ^^(2)^(2)^(1)^(1) T 
yy-l — A 12 -^23 -^31 +-^12-^23-^34-^41 -J^all 

w - F(2)\-(i)y(i) ^ Tr(i)3r(2)y(i)y(i) T 

^3 — -^12-^23-^31 +-^12-^23-^34-^41 ~ -l^ all 
W - F(2)\-(2)Tr(l) _L v(l)3r(l)y(2)y(l) T 

W4 — A;^2 -^23 -^31 +-^12-^23-^31-^41 ~ ^ all , {O.Z) 

where Taii x'^ x'^ x'^ X^^ x'^^ is composed of a full Eulerian cycle. These superpo- 
tentials together with the quiver diagram create all possible TCYQTs for this diagram, and so these 
theories encode all possible toric moduli spaces. Thus the program has produced a complete set of 
new physical theories for the above cases. Further complete sets of TCYQTs have also been produced 
for larger x and giving many more quiver theories than previous known. 




Figure 13: An example of {G,E) = (4,8) quiver diagram, for which we present all superpotentials in 
Equation (|5.3p . 
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As a final example, below are tables of the full sets of possible superpotentials for one of the sixteen 
(4, 8) diagrams above, drawn now in Figure [T51 



W — — — Xl2-'f21-^23-''^34-'^42 + ^l,3-'^21-'^34-^42 

W — ~Xi^3X^^X4i + Xi2X23Xl4Xii + Xl,3X2lX|4X42 — ^12-^^21 -'^23-^^34 -^42 

W — X12X23X34X41 — Xi2^1^3X2lX|4X4i — ^23-''^34-'^42 + -'^l ,3-^^21-'^ 34-''^42 

W — —X12X1, 3X21X34X41 + X12X23X34X41 + Xi^3X2lX|4X42 — X23X34X42 

W = Xi^^X^^X^i — X12X23XI4X41 — Xi 3X21X34X42 + X12X21X23X34X42 

W — Xi^3X|4X4i — Xi2Xi_3X2lX|4X4i + X12X21X23X34X42 — X23X34X42 

W — — Xi2Xi_3X2lX34X4i + Xi^3X|4X4i + X12X21X23X34X42 — X23X34X42 

W = — Xi^3X34X4i + X12X23XI4X41 + Xi^3X2iX34X42 — X12X21X23XI4X42 

W — — Xi2Xi^3X2iX34X4i + X12X23XI4X41 + Xi^3X2iX34X42 — X23X34X42 

W = — Xi2Xi^3X2lX34X4i + X12X23XI4X41 + Xi^3X2lX34X42 — X23XI4X42 

W = Xi^3X34X4i — Xi2Xi_3X2iX|4X4i — X23X34X42 + X12X21X23XI4X42 

W = Xi^3X34X4i — Xi2Xi^3X2iX34X4i + X12X21X23XI4X42 — X23XI4X42 . (5-3) 



We see that there are a total of 12 inequivalent possibilities. 
5.3 The Space of Quiver Theories 
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3/0 


8/47 


28/330 


72/1277 


208/3941 
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3/0 


12/0 


44/645 


172/4888 


6 














4/- 


16/- 


84/- 



Table 1: Number of Quiver Diagrams/Quiver Theories with non-vanishing Abelian superpotential, plotted 
against E, the number of edges and x — E — G, the number of superpotential terms. 

It is now expedient to summarise our results, which we plot in Table [1] containing the total 
number of quiver diagrams and TCYQTs for a given x ~ E — G and E. It is seen that the lower 
bound on the edges for the possible existence of TCYQTs is obeyed by the numbers in the table, 
as well as the condition that x > 3 for there to exist a TCYQT having a non-vanishing Abelian 
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superpotential. To further confirm the results the numbers oi x = 2 quiver diagrams are compared 
with the generating function found in [16], which gives the number of base quiver diagrams that have 
at least one non-vanishing toric superpotential. The expanded function is: 

g{x) = 1 + 2x + + g^3 ^ g^4 ^ ^3^5 ^ ^j^6 ^ 24a;^ + . . . , (5.4) 

with the number of {N + 2, N) quiver diagrams given by the coefficient of the x^^"^^ term. 

We see immediately that there seems to be a discrepancy with the (4,6), (5,7) terms and indeed 
the rest of the x — 2 quivers. This is because [16] does not restrict the fields to be complex valued 
and even for vanishing Abelian superpotential the gauge ranks could be promoted to > 1 and the 
two-term superpotential becomes non- vanishing. This means that all diagrams in Figures [7] and [8] 
except the ones with red nodes admit a superpotential as the order of the fields in a superpotential 
term matters. The quiver diagrams with red nodes do not because an overall trace is taken on the 
superpotential, and the two terms in the superpotential correspond to the same Eulerian path and 
so are cyclically related. Therefore, one more quiver diagram is produced for the (4, 6) quivers and 
two more for the (5, 7) quivers and the results are then in perfect agreement. The number of quiver 
diagrams found must then be an upper bound on the generating function coefficient and we see that 
this holds for the new data. 

It is worthwhile to consider the trends of the numbers in Table [1] Although there is limited data, 
there are a number of observations that can be made which provide evidence for the algorithm running 
correctly, and give more insight into the space of TCYQTs. The first is that large x and small E do 
not produce very many diagrams. This is fairly intuitive, as there are not many ways to place edges 
between a small number of nodes. There are many more TCYQTs than quiver diagrams for high 
E and x, but this is expected, as there should be many ways of decomposing a quiver diagram into 
loops, and consequently many term collections. 

One of the most interesting aspects is that whilst the number of diagrams increases monotonically 
with E and x, the number of TCYQTs does not for low E. It appears the number of quiver diagrams 
increases faster with E, and the number of quiver TCYQTs increases faster with x, but it is hard to 
draw concrete conclusions from such a small data sample. However, we can at the very least say that 
the number of quivers becomes large quickly with E, and gets larger with E faster for larger x. 

As well as the production of the quiver theories, one of the main results of this paper is the 
algorithm itself. Each routine was tested individualljl*]and that, together with the production of 
correct quiver diagrams and TCYQTs for low x and E provides sufficient evidence to be confident 
that the algorithm is robust. By a rough complexity analysis, given in Appendix [B] we find the 
algorithm has complexity ^ 0{E°''^^^) for some positive constants a and b, which for a given x is 
polynomial in E. Examining the timing table in Appendix [b1 we see that this complexity roughly fits 

*The testing strategy and more algorithmic details are detailed in Appendix IB] 
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the data. Therefore it is feasible to use the program to produce quiver diagrams in a short time if the 
number of contractions is smaU, say a; < 7, and for large x if time is not an issue. 

6 Conclusions and Prospects 

We have produced complete sets of toric quiver theories for low orders and a method to calculate all 
possible TCYQTs for higher orders. Many of the generated quivers have appeared in the literature 
for the first time, providing us with large number of new theories to study. 

The idea of collapsing polygons to produce quivers and AB notation are new concepts. Moreover, 
the order parametre x = E — G, the number of terms in the superpotential, is also the number of 
contractions in this language, have allowed us to find a constructive solution to the x = 2 problem, 
and, in general, an efhcient algorithmic solution for arbitrary x. The AB notation is more efficient and 
more descriptive than the incidence matrix, in that it contains less redundant information, records 
the placement of self-adjoint loops and reflects the symmetries of the quiver diagram. 

We have shown that the toric condition has a combinatoric interpretation of splitting the graph 
into disjoint loops, which we have exploited to generate superpotentials by constructing all sets of 
disjoint loops. It is easy to see from this that all toric spaces are Calabi-Yau spaces, as all Eulerian 
graphs can be split into disjoint loops, the proof of which we again leave to Appendix El The solution 
suggests that it is a better idea to characterise the order of quivers not by their edges and nodes [E, G) 
as was in [16], but, rather, by the number of edges and the number of contractions {E, x), because it 
is X not and G which determines the number of elements in the AB. 

The computational complexity of the algorithm is roughly 0(iJ''^+^) for some positive constants 
a and 6, which is much better than the brute force case. It is unlikely that there exists an algorithm 
which has purely polynomial complexity in both E and x, though it may be possible to improve the 
algorithm by using one of the more complicated graph isomorphism tests available. Our plenitude 
of data can be used to study the space of general toric quiver theories more extensively than what 
has been possible before. As the algorithm is exhaustive, this data can be analysed systematically. 
Producing all these types of quiver theories would involve expanding the algorithm to take matrix 
valued fields into account when producing superpotentials. 

The toric quiver theories have been produced specifically so the 'forward algorithm' [32] can be 
applied in order to scan through lower order theories, giving us a new method to study the duality 
manifested in the AdS/CFT correspondance. In particular, we should go through all the quivers 
together with their superpotentials and start listing the VMS by explicitly solving for the D-term and 
F-term equations. The result should be all possible, stepwise, toric Calabi-Yau threefolds. Likewise, 
we could add Chern-Simons labels to the nodes and follow the extended forward algorithm of [16]. 
The results should be all toric Calabi-Yau fourfolds. This is of tremendous interest because so far we 
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do not have an inverse geometrical engineering method of constructing M2-brane world-volume gauge 
theories and exploring this space of theories under an algorithmic light should provide new insights 
and dualities. 
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Appendices 

A Some Technical Results on Quiver Constructions 

In this appendix we collect some technical results which are used in the paper. 
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A.l Shortest Distance and the Adjacency Matrix 

There is a relation betwen shortest distances between two nodes and the adjacency matrix of the 
graph. Consider the n*'' power of the adjacency matrix A = aj. Using summantion notation, We see 
that A'^ = a^aj. and so 

A- = a^^4l..4:_^al = {air . (A.l) 

So, since the adjacency matrix is one if there is an edge between i and j, and zero otherwise. Then 
(a^)" is only non-zero if there is some path that goes from i to j in n steps. As we require all the 
al''+'' to be non-zero, simultaneously, in the sum. So there must be a path from i to j of n steps in 
length if this element of A" ^0. The first time this is non-zero will be the shortest path between i 
and j, as this multiplication tries all routes. 

A. 2 A Quiver Admits an Eulerian Cycle 

As stated in the main text, every toric quiver admits an Eulerian cycle. Here we present the proof, 
due originally to Fleury in [24]. 

Proof. Consider a {E, G) quiver diagram, where E — G = x. Pick any node in the diagram. Now 
travel from this node, through the graph, never using any edge more than once. Every node we arrive 
at, we can then leave, apart from the node we started at. So we must end up back at the start, and 
have created a loop in the graph. Now remove that loop from the graph. The remaining graph is still 
a quiver, and we can perform the same operation. Continuing in this manner we can include every 
edge in some loop in the graph. Now reconstruct the whole graph, remembering the loops. Starting 
from the intial node, if we come to a node where another loop begins, we travel around it. In this 
way we travel through every edge once, and form an Eulerian cycle on the quiver. □ 

A. 3 CycUc Permutation of A and B 

There is a permutation symmetry on the A and B maps on the quiver. Consider an [E, x) quiver 
A = {ai...ax},B = {bi...bx}, and rotating this quiver to the next node which a contraction starts 
at. The first contraction is then 02, and then next step, 62, and so on until we reach a^- Then the 
previous last step b^, steps to the first contraction oi and 61 steps to 02. So a rotation of the quiver 
simply amounts to an identical cyclic permutation applied to both the A and B. 

A. 4 Mirror Symmetry of AB 

The mirror symmetry of a quiver is also reasonably nicely manifested. To construct the mirror image 
of a quiver, and remembering that we always start on a contracted node, we first contract ai. Then 
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we travel + ai steps to the next contraction. We then contract E — ax, and travel bx-i + a,x to the 
next contraction, then contract E — ax-i and so on. This gives a nice condition on the a's. That the 
sum of the a's is less than Ex/2, if we allow the Vs, to take any value. Because of this allowance on 
the 6's the AB can 'wrap around' the quiver more than once which is why we cannot use this and the 
moduli rules together. 

A. 5 Stopping Redundancy 

Setting 6 > 1 stops a redundancy when coupled with the moduli rules as in order to produce a 
redundant a then during the walk to collapse the polygon, we must arrive at a loop already formed, 
and then apply C{a) where a is the number of edges in the loop. If 6 > 1 we can never arrive at such 
a loop, as we can only create loops where we are, and then if we move forward once we cannot go back 
to the beginning of the loop. However, this doesn't stop us from making any possible quiver. As this 
c-node is now in front of us too, so we can contract any node in front of us with it, except for itself. 
In other words, this condition removes the redundant possibility whilst keeping all the ones that are 
not. 

B The Algorithm 

This appendix contains a brief set of documentation explaining how to test the algorithm. It also 
contains details of the implementation of the algorithm in C. The list below splits the program into 
its three sections, detailed in the flowchart in Figure [5l It contains the information on what each of 
the implemented routines do, and how they were tested. Following this, results are shown for the 
runtime of the algorithm and the calculation of the approximate complexity algorithm is also shown. 

B. l Implementation 

The isomorphism test works by creating a 'path' object which walks through the quiver. First we 
generate a quiver diagram from the AB we are testing. We walk all possible Eulerian cycles of the 
quiver to find all possible AB. The walker starts at a c-node in the diagram, and then starts walking, 
every time the walker has a choice of where to go, it recognises it is at a c-node and so a contraction 
must begin from this node, it starts then counting from this node until it returns to it, at which point 
it has a value for that contraction. After leaving a c-node, the walker counts how many edges it walks 
until it reaches the next c-node. This is a B value. The program walks all paths simultaneously and 
so, when it reaches a c-node, where there is a choice, the walker copies itself, and walks both paths. 
This creates a web of Eulerian cycles across the quiver, each with a certain AB value. We then pick 
the smallest AB of these and compare it to the AB we used to generate the quiver. If it is the same, 
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then we keep the AB, if it is smaller, we discard our tested AB. In this way, we find a representative 
of each graph isomorphism class and so generate all the quiver diagrams we require. 

To output the data, a routine for converting the AB into an incidence matrix is needed. An 
incidence matrix for the original polygon is created, then a list of contracted is made. This array 
is then used to sum the correct rows from the polygon incidence matrix, equivalent to contracting the 
set of nodes stored in the row. Wo can then form the incidence matrix by taking each of these rows. 
The incidence matrix is the ouput to a file, in a format suitable for the toric analysis of the graphs. 

To produce all the quiver diagrams, without the superpotential generation, a stripped down version 
of the program is used. This version does not contain the superpotential generation routine and has 
a very slightly different structure. As the supcrpotcntials arc not calculated, the self-adjoint loops 
become irrelevant, as they can go anywhere, and do not add any useful data to the diagram. So, a 
routine which checks for self-adjoint loops is used, which decreases the number of AB which have to 
be checked for isomorphism. This check occurs when the program creates an adjacency list, as it is 
easy to see if the loops are formed, as a node will be contracted with one next to it in the array. 

To generate the stripped down diagrams, the ones missing self-adjoint loops, the algorithm is called 
again, but for {N, E — 1), then for {N, E — 2), until E = N. This has the added benefit of producing 
all the quivers with the same N down to {N, N) for free. 

B.2 Break of Sub-routines in Quiver Generation 

The routines break down as follows: 
1. Counting section: 

(a) Main - Accepts intial variables and sets up the problem, tested by printing the input 
variables to the screen, and making sure that the routine calls the following ones. 

(b) Generate Quivers - Sets initial A and starts counting, tested by printing the intial A and 
checking that the final A is the first one which ovcrfiows. 

(c) Aperm - Counts through A's, tested by printing each A to the screen and checking that 
they are being counted through correctly. 

(d) Bperm - Counts through B's. tested by printing each B to the screen and checking that they 
arc being counted through correctly. Due to the increased complexity of counting through 
the B, due to the varying end points for the different b's, each ovcrfiow was checked by 
hand for the small E and x. Then, a testing routing was used to check that the AB never 
broke the moduli rules for the other higher order AB. 

(e) Incidencegcn - Calls the isomorphism checking routine for the quiver, produces incidence 
matrix, and outputs to file, calls superpotential generation routine. Tested by checking the 
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correct incidence matrix was produced for low complexity quivers, and it's easy to look in 
the file to check that the output is correct. 

(f) aadder - Handling routine, increments A, tested by printing the A before and after incre- 
mentation. 

(g) badder - Handling routine, increments the first b of B. Tested by printing the B before and 
after. 

(h) recurseadd - Handling routine, increments B, tested by printing the B to the screen. 
Isomorphism check section 

(a) Finalcheck - Creates all possible paths, steps the paths through the graph by calling the 
stepping routing and calls the comparison routine. Tested by printing all the created paths 
to the screen, by printing that a path has been stepped to the screen and printing the result 
of the comparison. 

(b) adjmaker - Creates the adjacency list for the given AB, tested by printing the adjacency 
list to the screen and comparing it to the corresponding graph in AB notation. 

(c) disconnectcheck - Checks to see if a path has made the graph disconnected, tested by 
printing the path, and the result of the disconnect check, and making sure that it was only 
giving the positive result for actually disconnected paths. 

(d) pathstepper - Steps a path along to the next adjacent nodes, creates more paths if the node 
is multiply connected, tested by printing the path to the screen every step, and checking 
that the steps are working correctly. 

(e) pathcompare - Compares all produced paths to the given AB, tested by printing the whole 
list of paths, to check their AB and comparing them by hand to the given AB. 

(f) deletepath - Deletes a path from the pathlist, tested by printing the list before and after 
the path has been deleted. 

(g) duplicatepath - Duplicates a path and puts it in the list before the path that was dupli- 
cated,tested by printing the list of paths before and after duplication. 

(h) initialisepath - Initialises all the variables in an empty path, tested by printing the path 
after initialisation. 

(i) quicksort - Handling routing - Sorts a list of data, tested by printing the list to be sorted 
before and after. 

Superpotential Generation Section 

(a) Constructsuperpotentials - Calls gcnerateterms, and constructs supcrpotcntials from the 
produced term collections. Tested by printing all term collections and putting them to- 
gether by hand to check they match up with the ones the program produced. 
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(b) cancclchcck - Checks to sec if two term collections contain any like terms, tested by printing 
the list and making sure all term collections are different. 

(c) printsuper - Prints the superpotential to a file, tested by reading the file. 

(d) assignmentadder - Changes the edge assignments in a systematic way, tested by checking 
that all assignments are tested, by working out what they should be by hand, and printing 
the program generated ones to the screen. 

(e) generateterms - Generates all the terms for a given edge assignment and adds them to the 
list of term collections, tested by taking some edge assignments, and running through the 
loops by hand, then comparing them to the term collections printed on the screen. 

(f ) identcheck - Checks whether any two term collections are the same, tested in the same way 
as cancelcheck. 

(g) deletecol - Deletes a term collection from the list, tested by checking the pointer to the 
term collection is made null after the term collection is deleted. 

(h) termsort - Takes a term collection and organises it so the numbers of fields, and the entries 
in a term are in order, tested by printing the terms before and after and checking that they 
are put in order. 

(i) tquick - Does the sorting of the terms, tested in the course of testing termsort. 
(j) fquick - Does the sorting of the fields, tested in the course of testing termsort. 

Each routine was tested to check that it produces the output that it is supposed to, and routines 

were tested in such a way that any dependencies were tested before the routine which was dependant 
on them. Testing the counting section was a simple matter, as that had to be done was to check 
that the program was indeed counting through all of the potential quivers, which we could look at 
the screen and see. The isomorphism check section was tested by dumping the contents of given 
path structures to the screen, so we could check that the adjacency list and pathstcppcr were working 
properly. The disconncctchcck was tested by following paths which became disconnected and checking 
that the checker picked them up. Superpotential generation was tested by dumping the contents of 
given term collections to the screen and counting the number of edge assignments that should be 
counted through. 

Obviously the ideas behind the algorithm work, and we have proven this earlier on in the main 
text. As a test, we applied the algorithm, by hand, to a number of sets of low order quivers, and 
compared the results with the ones the program gave. The results it gave were indeed the same. 
However, it is hard to probe any further as it becomes incredibly laborious, and a person is more 
likely to make mistakes than a computer is at that stage. The algorithm itself however is quite robust 
and quite modular. 
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We were also able to compare the results of the program in the x — 2 case to the direct implemen- 
tation of the solution to the warm-up problem and check that both programs produce the same results, 
which they do. To give a better idea of exactly how some of the routines work, the isomorphism test 
and output routine are covered in more detail below. 

B.3 Computational Results 

To analyse the efficiency of the program, we look at the below table, which records the time taken 
for the program to run. We see that the program is "fast" in E, and "slow" in x. This fits with our 
analysis of the complexity of the algorithm below. 

Table 1: Table showing time (in seconds) for Quiver Diagrams/Quiver Theories generation 
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v- 
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2/- 


1/1 








6 


2/- 


2/5 


4/4 






7 


2/- 


2/5 


2/7 


11/34 




8 


3/- 


3/6 


7/25 


9/151 


15/- 


9 


3/- 


8/15 


27/117 


63/908 


16/- 


10 


4/- 


24/45 


147/571 


835/6690 


1895/- 



Simple rules to note are that loops are polynomial in the number of entries in the loop, and that 
counting through numbers is exponential in the number of digits. This is a very rough analysis and 
should just give an idea of how efficient this algorithm is, i.e., whether it is fast or slow. 

We see that the counting section of the program is exponential in its complexity in a:. As it counts 
through quivers the order of magnitude of the number of quivers it has to count through is roughly 
0{E^). There is, as far as we can see, no way around this as every quiver must be checked. The 
isomorphism checker has complexity 0(i?2'^y, as has already been discussed. The superpotential 
generator has to count through all assignments and generate the terms for each one. Term generation 
is polynomial time in E as it simply loops through all the edges and records the term collections. 
However, the assignments are again exponential in x, so we get a complexity like 0{E°''''^^). 



*As discussed in [26]. This seemingly amounts to a brute force method, however, because we only start paths on 
the c-nodes the method does in fact scale polynomially with E, as we have 2^ paths maximum, which we step E times. 
Hence the number of steps is E2^ 
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